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Single-photons play an important role in emerging quantum technologies and information process- 
ing. An efficient generation technique consists in preparing such states via a conditional measure- 
ment on photon-number correlated beams: the detection of a single-photon on one of the beam can 
herald the generation of a single-photon state on the other one. Such scheme strongly depends on 
the heralding detector properties, such as its quantum efficiency, noise or photon-number resolution 
ability. These parameters affect the preparation rate and the fidelity of the generated state. After 
reviewing the theoretical description of optical detectors and conditional measurements, and how 
both are here connected, we evaluate the effects of these properties and compare two kind of devices, 
a conventional on/off detector and a two-channel detector with photon-number resolution ability. 



I. INTRODUCTION 

A driving component for the development of diverse 
quantum information applications is the ability to effi- 
ciently engineer non-classical states of light In par- 
ticular, a great number of communication or computing 
protocols require single-photon states [2j . The control of 
single quantum objects provides a way to generate such 
states @]. Another procedure, probabilistic but heralded, 
consists in measuring one mode of an entangled state, 
which results in projecting the other mode according to 
the measurement result: in the ideal case of twin beams, 
the detection of a single-photon in one beam projects 
the other one into a single-photon state The physi- 
cal realization of this scheme critically relies on the per- 
formances of the twin-photon generator but also of the 
single-photon detector used to herald the preparation. 

Indeed, realistic detectors do not perfectly resolve pho- 
ton numbers, have limited quantum efficiency and can be 
noisy [5|. These properties directly translate into the fi- 
delity of the generated state and its preparation rate. 
The goal of this paper is to evaluate how these proper- 
ties impact the conditional preparation of single-photon 
states. We will consider the typical case of twin beams 
produced by two-mode spontaneous parametric down 
conversion. To illustrate the effect of the heralding de- 
tector properties, we will analyze two different detectors, 
a conventional APD and a single-loop time-multiplexed 
detector. The present work provides a detailed analysis 
of the experimental characterization of these detectors @ 
and of the modeling of the single-photon state prepara- 
tion using such detectors. 

The paper is organized as follows. In section [TTl we 
first recall the general description of optical detectors by 
positive operator valued measure (POVM) and describe 
the conditional measurement scheme based on entangled 
beams. Given the detector POVM, we derive the gen- 
eral expression of the generated state. We illustrate this 
section by considering two specific detectors. Section HTT1 



is then devoted to the experimental characterization of 
such detectors via quantum detector tomography. In Sec- 
tion |lVl by numerically simulating the entanglement re- 
source, we predict the properties of the generated state 
that would be obtained if such detectors are used. We 
compare the two devices. We finally discuss the depen- 
dence of the generated state in the general case obtained 
by using the theoretical expressions of the POVM ele- 
ments. This study enables us to determine the fidelity 
and preparation rate for the two heralding detectors, as 
a function of the noise background, quantum efficiency 
and entanglement resource. 



II. THEORETICAL DESCRIPTION OF 
SINGLE-PHOTON DETECTORS AND 
QUANTUM STATE ENGINEERING 

Quantum state engineering is an interplay between the 
entangled resource and the measurement performed by 
the heralding detector. In this section, we remind how 
the detector properties can be described in the most gen- 
eral case by a positive-operator valued measure and we 
provide the general expression of the state prepared by 
a measurement performed on one beam of the two cor- 
related ones generated by spontaneous parametric dow- 
conversion. Two particular detectors are considered. 



A. Description of a detector by a POVM 

A detector can be fully described by its positive op- 
erator valued measure pj, i.e. by a set of semipositive 
hermitian {n„} associated to each of its N possible out- 
comes labeled by integer "n" . These operators sum up 
to the identity, J2n=o ^ n = They enable to make 
predictions on the measurement outcomes. For any im- 
pinging states p on the detector, the probability p n (p) of 
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obtaining the outcome "n" is given by: 
Pnifl) = Tr[/5n„] 



where Tr stands for the trace. 

In the following, we will consider the case of photon 
counters, which are phase-insensitive. In this specific 
case, the outcomes "n" correspond to "n clicks". For 
ideal detectors with perfect photon number resolution 
ability, the number of clicks directly gives the number of 
impinging photons. Therefore, the corresponding POVM 
is the set of projectors on the photon-number states: 
II„ = \n)(n\. 

However, realistic detectors are affected by non-unit 
quantum efficiency, imperfect counting ability and in- 
trinsic noise. The general expression for Tl n can then 
be written as a sum of projectors: 



n„ = y"Vfc, n |fc) (k\. 



k=0 



Let us note that the coefficient r^, n can be read as: 



As stated by equation (J3j) , the coefficient rk, a ff is given 
by the expectation value of the operator Yl ff for an im- 
(1) pinging Fock state |fc). In the absence of noise, this value 
is equal to the probability of not detecting k photons for 
a quantum efficiency rj: it is indeed equal to the prob- 
ability to not detect a photon, 1 — rj, to the power of 
k. 

The second detector we consider is a single-loop 
time-multiplexed detector with photon-number resolu- 
tion ability (Fig. [TJ:). Indeed, it is able to provide an 
additional information by distributing the light to be de- 
tected in two temporal bins. The incoming photons are 
sent to a beamsplitter (reflectivity R in intensity) whose 
outputs are delayed by a time interval AT and recom- 
bined before being detected by an APD Q. If two pho- 
tons impinge on the detector, they are split toward dif- 
ferent paths with probability 2R(1 — R) and, in this case, 
they are detected by the APD as two subsequent events. 
Three outcomes are thus possible: no click, 1 click and 
2 clicks. For a given quantum efficiency rj of the APD 
used for the detection of the two bins, the POVM ele- 
ments of the TMD take the general form: 



(2) 



>'k. 



Tr[n„ \k)(k\] = P{n\k). 



(3) 



It thus corresponds to the conditional probability P(n\k) 
to obtain the outcome "n" giving an impinging Fock state 
\k). This expression is very convenient to determine the 
POVM coefficients of a given detector implementation. 



B. General expression of the POVM for two 
experimental cases 

We will consider here the case of two single-photon 
counters: a conventional on/off avalanche photodiode 
(APD) and a home-made single-loop time-multiplexed 
detector (TMD). We give the expression of the POVM 
elements in the general case, including non-unit quan- 
tum efficiency rj and dark noise with Poissonian statistics, 
which corresponds to many practical situations. 



1. Including non-unit quantum efficiency 

A conventional APD has a limited counting ability as it 
provides a binary outcome: it is only able to tell whether 
photons (whatever their number) have been detected (on, 
one click) or not (off, no click). In the ideal case, the two 
corresponding POVM elements are given by Il ff = |0) (0| 
and n o „ = 1 — |0)(0|. In the case of non-unit quantum 
efficiency n, the POVM for the APD detector can be 
written as 



n off (v) = J2Zo^off(v)\k)(k\ = £r=o(l - V) k \k)(k\ 

H-oniv) = J2T=0 r k,on(v)\ k )( k \ = 1- fW(?7, V). 



H-o(v) = 12^=0 r k,o(v)\k)(k\ 
Mv)=T l k= rk,i(v)\k)(k\ 

[ n 2 ( ?7 ) = i-no(r ? )-ii 1 (r ? ). 



(5) 



The POVM coefficients are calculated following the def- 
inition given in ([3]): 

r fc , (r?) = P(0\k) = £ JZ"(1 - R) k -pQr k ^ off r p , off 



fe=0 



= (!-»?)* 



and 



rk 



xfo) = P(l\k) = j2B?(l-R) 



k—p 



(6) 



(7) 



X- (j"k—p,on fp.off H~ Tp.on ^k—p,off^) 



= (i-v) k pR p (i~R) k - p Q 



x(-2 + (1-77)*-* + (1-77)-") 

-) fc + (i + T^ 

rj' v 1 — rj' 



d-v) k [-2 + (i + -^) k + (i + ^L) k 

I x 1 — ri' v 1 — n' 



(4) 



2. Including dark noise 

We now include the effect of Poissonian dark counts 
in the description. We note v the mean number of dark 
counts in the detection window, which is defined before 
the 'black-box' detector: e~ v is thus the probability of 
no dark counts whatever the specific detector implemen- 
tation. For instance for the TMD, where the mode is 
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split into two temporal bins, e~ v l 2 is the probability of 
no dark counts into each bin and [e~ v l 2 ) 2 = e~ v for the 
two bins. 

The general expression of the POVM in the presence of 
noise can be written as a function of the POVM without 
noise. Using <j3j) , one obtain for the APD 



rk,on(v,v) = 1 - e~"r k>off (rj) 



while for the TMD it can be read as: 



(8) 



rk,o{r),v) = e u r k ,o(r)) rq s 
r M (r?, v) = e- v l 2 r k>1 {r,) + 2{e^l 2 - e-")r fc ,ofa). { 1 



C. Conditional state preparation with twin beams 

We now detail the conditional preparation. As ex- 
plained before, a general strategy for the conditional 
preparation of quantum states relies on the use of en- 
tanglement: due to quantum correlations, measuring one 
mode of an entangled state results in projecting the other 
mode according to the result of this measurement. By 
denoting pab the two-mode state, a measurement per- 
formed on mode B characterized by the POVM IT n leads 
to the preparation on mode A of the conditional state p c , 
which can be written as: 



1 



Pc 



Pr(n) 



Tr, 



pab i a ® n r , 



(10) 



where Tr^ stands for the partial trace over mode B and 
i-A is the identity operator acting on mode A. Pr(n) is 
the probability to obtain the outcome "n" on mode B 
and is thus proportional to the preparation rate. It can 
be read as: 



Pr(n) = Tt ab p AB U ® H 



(11) 



In this paper, we consider as an entangled state re- 
source the twin beams produced by perfect two-mode 
spontaneous parametric down-conversion (SPDC): 



oo 



k=0 



\k) A \k)i 



(12) 



where the squeezing parameter A, which varies between 
and 1, is related to the intensity gain in the down- 
conversion process. This parameter determines the av- 
erage photon number per mode as A" = |A| 2 /(1 — |A| 2 ). 
From the previous general expressions (fT0|) and (fTTj) . it 
follows that the conditional state p c is given by 



Y.T^rr k ,n\k){k\ 
Y,T=o \M 2k rk,n 



(13) 



where r k n are the coefficients of the POVM element il„. 
The probability Pr(n) can also be written as a function 



of the squeezing parameter and POVM coefficients as 



Pr(n) = (l-|A| 2 )£|A 



2 k 



(14) 



fc=0 



For an ideal photon-number resolving detector, 
r k,n( r ]y 1 ') = 4,m the conditional state resulting from a 
heralding measurement II n = |n)(n| reduces thus to a 
Fock state |n) for any values of A as expected. Only the 
preparation rate varies with this parameter. In the more 
general case, p c tends to the normalized POVM element 
in the limit of high gain: 



A-H 



n„ 



Tr II 



(15) 



Note that defining the non-classical character of a 
quantum measurement is not straightforward as the 
POVM elements are not normalized and thus do not 
correspond to a quantum state. However, conditional 
preparation provides a good insight into this subtle ques- 
tion. Given equation (|15p. the non-classical properties 
of a detector can be seen as the non-classical proper- 
ties of the state generated by a conditional measurement 
in the limit of high gain. Significantly, this expression 
also corresponds to the expression of the state retrod- 
icted from the measurement outcome [Tol - fl2l | . This state 
enables to make predictions (called retrodictions) about 
state preparation leading to this result. The non-classical 
properties of the detector can then be associated in this 
approach with the properties of the retrodicted state [l2j . 



III. QUANTUM DETECTOR TOMOGRAPHY 
AND EXPERIMENTAL DETERMINATION OF 
THE POVM 

We now turn to the experimental characterization of 
detectors. The complete characterization of a detector 
requires, for each of its POVM elements ±l ra , the knowl- 
edge of all the coefficients {rk, n (r], v)}. This knowledge 
can be obtained experimentally by performing the so- 
called quantum detector tomography (QDT) [q. Il3l4l9| . 
We present here the experimental setup for QDT and 
the results for the two detectors we are focusing on in 
the present study. 



A. Procedure and setup 

The procedure for QDT consists in probing the de- 
tector by sending a tomographically complete set of 
known states {pj} and in recording the number of times 



(fn 



each outcome 



occurred in the measurement 



for each probe state. By using equation (fTJ , a maximum 
likelihood algorithm (ML) allows for reconstructing the 
POVM elements {II„} by searching for the set of the 
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AT = 150ns 



FIG. 1: Experimental quantum detector tomography, (a) The detector, shown here as a black-box, is probed with coherent 
pulses coupled into a single-mode fiber. A fraction of the light is tapped and measured with a calibrated power-meter to adjust 
the amplitude of the probe. Two detectors are characterized: (b) a single-photon counting module (APD) and (c) a two-channel 
time-multiplexed detector. 



fk,n{ r li 1 ') coefficients that provide the optimal matching 
between the measured f n .j and the theoretical p n ,j- Sig- 
nificantly, the ML does not require any a priori assump- 
tions on the detector under investigation. The details of 
the procedure can be found in 



In our experiment, we used as probe states a set of co- 
herent states |aj)(o!j| of different amplitudes, generated 
by strongly attenuating the light coming from a pulsed 
laser source. The experimental setup is shown in Figure 
1. Light pulses at 795 nm are provided by a femtosec- 
ond laser source (MIRA900) with a repetition rate set to 
1.187 MHz thanks to a pulse picker device (9200 Coher- 
ent). In order to tune the average photon number per 
pulse, \oij\ 2 , the laser beam is sent through a variable 
attenuation device consisting of a half wave plate (A/2) 
and a polarizing beam splitter (PBS). The power of the 
output beam is controlled by tapping a known fraction 
toward a power meter. A set of well calibrated optical 
densities (total O.D.&7) supplies the attenuation neces- 
sary to reach the regime of few photons per pulse: |aj| 2 
can be indeed adjusted from to 100 by steps of 0.1. The 
probe beam is sent to the "black-box" containing the un- 
known detector through an optical fiber. The losses due 
to the coupling are taken into account when estimating 
the average photon number. The error on the average 
photon number is of the order of 5%, mostly due to the 
residual laser intensity fluctuations between subsequent 
pulses. 



B. Experimental reconstruction in the absence of 
noise 

We performed the detector tomography procedure 
on the two different devices presented in section [II] 
a standard on/off avalanche photodiode (APD) and a 
single- loop time- multiplexed detector (TMD). The APD 
(Perkin Elmer SPCM-AQR-13) is sketched on figure lb 
while the multiplexed detector is shown in figure lc. For 
the TMD, the light is split toward two channels: pulses 
going to the upper channel undergo a delay of AT =150 
ns with respect to those in the other channel. The two 
paths are eventually spatially recombined using orthog- 
onal polarizations (this recombination avoids extra 50% 
losses) and sent to an APD. As AT is larger than the 
detector dead time, the arrival of two photons will cor- 
respond to two subsequent events. This setup leads to 
three possible responses: off, 1 click or 2 clicks. For both 
detectors, the overall quantum efficiency has been set to 
0.28 ±0.02. 

Figure 2 gives the experimental results of the QDT per- 
formed on both detectors (first line for the APD, second 
line for the TMD). Column (a) provides the frequency of 
the different outcomes as a function of the average photon 
number \a\ 2 . Column (b) then gives the reconstructed 
operator, IT on for the APD and 111 for the TMD, and 
column (c) finally shows the corresponding Wigner func- 
tions. As it can be seen, both are negative at the origin, 
which is a strong signature of the non-classicality of the 
detectors. Let us underline that this detector negativ- 
ity is a critical property for quantum state engineering. 
Indeed, the preparation of a quantum state with nega- 
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FIG. 2: Experimental data and reconstruction of the POVM elements. The first line corresponds to the APD while the second 
line gives the results for the TMD. Column (a) gives the bare tomographic data, i.e. the probability of the different outcomes 
as a function of the average photon number of the coherent probe pulses. Column (b) shows the reconstructed POVM II on for 
the APD and fii for the TMD. Finally, Column (c) provides the Wigner representation of these operators. 



tive Wigner function requires an heralding detector with 
negative Wigner function if the entangled resource has a 
positive one. 



0.30 

o 

S 0.10 
0.00 

0.15 

S o.io 

X 0.05 
0.00 
-0.05 









1 


•' v = 0.16 











|a| 2 



C. Experimental reconstruction in the presence of 

noise 

We performed similar QDT by adding noise to the de- 
tection process. The background noise is controlled by 
injecting on the detector, together with the probe states 
required for the QDT, a continuous- wave coherent beam 
at 1064 nm (Diabolo, Innolight). By adjusting the laser 
power, we are able to mimic different dark count levels 
with Poissonian statistics. As in section [Til we note v 
the mean number of dark counts in the detection win- 
dow. Figure 3 gives the Wigner functions for the two 
devices and for different noise levels. As it can be seen, 
the negativity decreases with the noise, and eventually 
disappears. This gradual disappearance of the negativ- 
ity witnesses in a quantitative way the degradation of the 
detector performances. 



IV. SINGLE-PHOTON GENERATION WITH 
APD AND TMD 



FIG. 3: Effect of dark counts. The cross-sections of the 
Wigner functions for the APD (a) and TMD (b) are given for 
different levels of dark counts. The inset of figure (a) gives 
the probability of the outcome 'on' for the different levels of 
noise while the inset of figure (b) provides the probability of 
the three possible outcomes (black='no click', red='l click', 
blue=2 click') for a given value of the noise background. 



We now turn to the conditional generation of single- 
photon states using these two specific detectors. By sim- 
ulating the entanglement resource, we first examine the 
state one would obtain using the characterized detectors. 
We then consider the more general case given by the theo- 
retical POVM coefficients derived in section|TTl in partic- 
ular for high quantum efficiency. This section illustrates 
the difference between the two detectors and the condi- 
tions in which the photon-number resolution ability of 
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FIG. 4: Effect of the dark counts on the conditional state 
generation. The state is calculated from Eq |13l by using the 
experimentally reconstructed POVM coefficients for the APD 
and for the TMD. Figure (a) gives the fidelity between a 
single-photon state and the state generated using the APD 
(dashed) or the TMD (straight), as a function of A and for 
different dark count levels. Figures (b) and (c) represent the 
cross-section of the Wigner functions of a single-photon state 
(green), of the states generated using an APD (red) and us- 
ing the TMD (blue), in the absence of noise (y = 0) and 
respectively for A = 0.1 and A = 0.7 



the two-channel TMD is useful. 



A. Fidelity of the state generated with the 
characterized APD and TMD 

Given the reconstructed POVM elements for the APD 
and the TMD, expression (|13[) enables us to determine 
the state one would conditionally generate by using 
these detectors and correlated photon pairs generated by 
SPDC. The fidelity of the generated state with the target 
one (pure state) can then be calculated as the overlap of 
the Wigner function, W c (x,y) associated to p r and the 
Wigner function of the target state W t (x,y) [20l l2lj|: 

p _ W c {x, y)W t (x, y)dxdy 
\W t \ 2 (x,y)dxdy 

In the following we calculate the fidelity with a single- 
photon state, whose Wigner function can be written as 
W t (x, y) = f e - 2 (* 2 +^[-l + 4(x 2 + y 2 )}. 

Figure 0^ gives the fidelity of the generated state for 
the APD and the TMD as a function of the squeezing pa- 
rameter A. We recall that increasing A leads to a higher 
number of photons in the twin modes. Therefore, in the 
limit of A — >• 1, the fidelity goes to due to multi-photon 
pairs, which can not be resolved by the present detec- 
tors. Very low squeezing parameter is usually required 
to overcome their limited photon-number resolution. In 
this case, the two detectors behaves similarly. Figures 
Hb and St provide the cross-sections of the Wigner func- 
tions in the absence of noise and for two different values, 
A = 0.1 and A = 0.7. In both cases, the generated states 
are very close. In the first case, however, the fidelity with 
a single-photon state is high while it strongly decreased 
for higher values. 
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FIG. 5: Fidelity and event rate in the absence of noise (u—0) 
and for a quantum efficiency ?7=0.9. Plot (a) gives the fi- 
delity between a single-photon state and the state generated 
using the APD (red) or the TMD (blue) as a function of the 
squeezing parameter A. Plot (b) shows the event rate R as a 
function of the targeted fidelity (A is chosen accordingly). 



In all the configurations, the TMD enables to reach a 
slightly larger fidelity due its photon-number resolution 
ability (yet limited). This gain appears for increasing 
value of A, while it is expected to disappear for vanishing 
A due to the the low-photon number occupancy of the 
modes. As it can be seen, the difference between the re- 
sults one would obtain using the APD and a single-loop 
TMD is however small. This is due to the low quantum 
efficiency of the avalanche photodiode used here. In the 
following, we present the general case, without being lim- 
ited by the efficiency 77 and noise v imposed by a specific 
experimental implementation. 



B. General case 

By using the theoretical expressions for the POVM co- 
efficients as given in section 2, we can predict the modi- 
fication of the conditional states and of their generation 
probability (the rate of heralded events) as functions of 
the detector parameters. 

In an ideal measurement, involving a projective de- 
tector the result 1-click of a measurement performed on 
mode B heralds the generation of a single-photon state 
in mode A (F =1) , whatever the statistics of the ini- 
tial twin beams, i.e. for all values of A. In this case, 
the squeezing parameter A only affects the probability of 
generating the state, Rideai = (1 — |^| 2 )|^| 2 - For real- 
istic detectors, the detector parameters affect both the 
generated state and the event rate. We compute here 
the fidelity of the generated state with a single-photon 
state, by making use of the expressions of rk, on and ri~ t i 
of equations (JHHJ). For the APD, the fidelity F on can be 
written as 

_ A 2 -(l- e -(l-?7)) (m 

b on - I e -„ (17) 

1-A 2 1-(1- J) )A 2 

and the rate of events R on is given by: 

e-(l-A 2 ) 
Ron ~ 1 ~ l-(l-r,)A2- (18) 
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FIG. 6: Fidelity in the presence of noise and for a quantum 
efficiency n=0.9. Plot (a) gives the fidelity between a single- 
photon state and the state generated using the APD (red) or 
the TMD (blue) as a function of the squeezing parameter A 
and noise v. Plot (b) represents the fidelity as a function of 
A for a noise i/=0.08. 



Similarly, for the state heralded by the TMD, we obtain 
a fidelity 

^ ^.(1-^-^) (19) 

l-(l-?)/2)A ;i l-(l-r))A 2 

with a corresponding event rate R\\ 

Rl = 2(1 ~ {l-(l- V )\2 + l-(l-r;/2)A 2 ) ' 

(20) 

As an example, we consider the case of a high quan- 
tum efficiency {q — 0.9). In the absence of noise (y = 0), 
figure [5] gives the fidelity and preparation rate as a func- 
tion of the squeezing parameter. As already seen in the 
previous section, the fidelity is decreasing with increas- 
ing A due to the multi-photon pairs emitted in the SPDC 
process. While the difference between the two detectors 
was small for a low quantum efficiency, it can be clearly 
seen now. The photon-number resolution ability of the 
TMD, yet limited, enables us to reach a larger fidelity for 
all values of A. Similarly, for a targeted fidelity, the TMD 
also provides a higher count rate, as shown in figure [SJd, 
as a higher A can be used. 

We finally compare the two detectors in the presence 
of noise. Figure [6] shows the fidelity as a function of 



A and v, still for a quantum efficiency 7/ = 0.9. For a 
given value of noise, the TMD enables to reach a better 
fidelity. As the noise increases, the difference between 
the two detectors decreases. As it has been underlined 
previously, for a given noise level, there exists an optimal 
value of A for the fidelity. For the squeezing parameter 
going to zero, the fidelity drops as the signal to noise 
ratio decreases, i.e. the noise results in a large rate of 
false heralded events relative to the low signal. 

V. CONCLUSION 



Beyond its fundamental significance, the study of 
optical detectors has become in the last years more 
and more crucial for quantum information, in partic- 
ular in the context of quantum state engineering and 
measurement-driven processes. We have presented here 
the general framework for conditional state preparation, 
which is an interplay between the resource and the 
heralding detector. In particular, we have taken the 
example of the conditional preparation of single-photon 
states and compared two detectors that can be used 
for such a preparation, an avalanche photodiode and a 
single-loop time-multiplexed detector. Four parameters 
are important in such a scheme: the entanglement 
resource defined here by the squeezing parameter A, i.e. 
the photon-number occupancy of the correlated modes, 
the photon-number resolution ability of the detector, 
i.e. its capability to discriminate multi-photon events, 
the noise of the detector and its quantum efficiency. We 
investigated the heralded state one would obtain given 
different sets of these parameters. For low quantum 
efficiency, as we shown experimentally, the TMD only 
provides a small gain in the heralding process. However, 
the effect of the photon-number resolution ability of 
the TMD clearly appears for larger quantum efficiency, 
enabling in particular a higher preparation rate for a 
given fidelity of the generated state. 
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